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Abstract Numerical simulations of magnctosonic wave formation driven by 
an expanding cylindrical piston are performed to get better physical insight 
into the initiation and evolution of large-scale coronal waves caused by coro- 
nal eruptions. Several very basic initial configurations are employed to analyze 
intrinsic characteristics of the MHD wave formation that do not depend on 
specific properties of the environment. It turns out that these simple initial 
configurations result in piston/wave morphologies and kinematics that repro- 
duce common characteristics of coronal waves. In the initial stage the wave and 
the expanding source- region cannot be clearly resolved, i.e. a certain time is 
needed before the wave detaches from the piston. Thereafter, it continues to 
travel as a so-called "simple wave" . During the acceleration stage of the source- 
region inflation, the wave is driven by the piston expansion, so its amplitude 
and phase-speed increase, whereas the wavefront profile steepens. At a given 
point, a discontinuity forms in the wavefront profile, i.e. the leading edge of the 
wave becomes shocked. The time/distance required for the shock formation is 
shorter for a more impulsive source-region expansion. After the piston stops, the 
wave amplitude and phase-speed start decreasing. During the expansion, most 
of the source region becomes strongly rarified, which reproduces the coronal 
dimming left behind the eruption. On the other hand, the density increases at 
the source-region boundary, and stays enhanced even after the expansion stops, 
which might explain stationary brightenings that are sometimes observed at the 
edges of the erupted coronal structure. In addition, in the rear of the wave a 
weak density depletion develops, trailing the wave, which is sometimes observed 
as weak transient coronal dimming. Finally, we find a well defined relationship 
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between the impulsiveness of the source-region expansion and the wave ampli- 
tude and phase speed. The results for the cylindrical piston are also compared 
with the outcome for a planar wave that is formed by a one-dimensional piston, 
to find out how different geometries affect the evolution of the wave. 

Keywords: Waves, Magnctohydrodynamic; Waves, Shock; Corona; Coronal 
Mass Ejections; Flares 

1. Introduction 

Explosive expansion of coronal structures associated with CME/flare eruptions 
frequently creates large-scale large-amplitude waves and shocks in the solar 
corona (for recent reviews, presenting detailed overview of various aspects of 
this phenomenon, see Warmuth, 2007; Vrsnak and Clivcr, 2008; Wills-Davey 
and Attrill, 2009, Warmuth, 20f0; Gallagher and Long, 2011; Zhukov, 2011; 
Patsourakos and Vourlidas, 2012). These global disturbances are observed as 
EUV coronal waves, chromospheric Moreton waves, type II radio bursts, moving 
soft X-ray, and/or radio sources (see Warmuth et al., 2004a; Vrsnak et at, 2006; 
Olmedo et al., 2012), as well as sharp fronts in the white-light coronagraphic 
CME images (e.g. Ontiveros and Vourlidas, 2009). In recent years, this phe- 
nomenon was the subject of many studies that focused on various observational 
and theoretical aspects, including the morphology, kinematics, source-region 
characteristics, shock formation, three-dimensional propagation, etc. (for a brief 
overview of recent research activities see, e.g., Section 8 in Klimchuk et al, 2009 
and Section 9 in van Driel-Gesztelyi et al., 2012). 

At low coronal heights, where disturbances are observed in the EUV range, 
waves usually become recognizable at a distance of rj 100 -200 Mm from the 
source active region (e.g. Vcronig, Temmer, and Vrsnak, 2008; Patsourakos et al., 
2009; Kienreich et al, 2011; Muhr et al, 2011). Thus, EUV waves are observed 
while propagating through a quiet corona, where the magnetic field is predom- 
inantly vertical. Consequently, a low-coronal wave segment can be considered 
as a perpendicular magnctohydrodynamical (MHD) wave (magnetosonic wave). 
Typical velocities of EUV waves are a few hundred kms" 1 (for details see, e.g., 
Thompson and Myers, 2009; Warmuth and Mann, 2011 and references therein). 

New detailed observations reveal that the wave amplitude initially increases 
and at the same time the wave accelerates. Eventually, after a phase of ap- 
proximately constant speed, the wave decelerates to velocities typically around 
200-300 kms" 1 (e.g. Long et at, 2008; Muhr et at, 2013; Temmer, Vrsnak, 
and Veronig, 2012), where faster waves show a stronger deceleration (Liu et al., 
2010; Kozarev et al, 2011; Ma et al, 2011; Warmuth and Mann, 2011; Cheng 
et al., 2012; Olmedo et al., 2012). It was also found that waves of higher speed 
have larger amplitude (Kienreich et al., 2011). During the constant-speed and de- 
celeration stage, the amplitude of the perturbation decreases whereas its profile 
broadens. Such behavior is usually interpreted as a typical signature of a freely 
propagating "simple wave" (for terminology we refer to Vrsnak, 2005; Warmuth, 
2007). 
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The fastest waves are frequently accompanied by type II radio bursts (e.g. 
Klassen et al, 2000; Bicscckcr et al, 2002; Warmuth et at, 2004b; Veronig et ai, 
2006; Vrsnak et ai, 2006; Muhr et ai, 2010; Ma et al, 2011; Kozarev et al, 
2011), which reveal the formation of a coronal MHD shock. Such waves may 
also generate Moreton waves (e.g. Warmuth et al, 2004a; Vrsnak et ai, 2006; 
Muhr et al., 2010; Asai et at, 2012; Shcn and Liu, 2012) if the pressure jump 
at the shock front is strong enough to push the inert chromosphcric plasma 
downwards, i.e. if the shock amplitude is high enough. 

Generally, coronal waves and shocks could be generated by the source-region 
expansion either related to a coronal mass ejection (CME), or a pressure pulse 
caused by the flare-energy release (for a discussion see Vrsnak and Oliver, 2008). 
Whereas in many events the source-region expansion could be clearly identified 
with the impulsive-acceleration stage of a CME (e.g. Patsourakos, Vourlidas, 
and Stcnborg, 2010; Veronig et ai, 2010; Grechnev et ai, 2011; Kozarev et ai, 
2011); in some cases there are indications that the shock is initiated by a flare 
(e.g. Vrsnak et ai, 2006; Magdalcnic et ai, 2010; Magdalenic et ai, 2012). 

Whatsoever the driver is, perpendicular MHD shocks are created by plasma 
motion perpendicular to the magnetic field. For example, a supersonic motion 
of small-scale ejecta would produce a shock (see, e.g., Klein et al, 1999), in a 
similar manner as that in which supersonic projectiles create shocks in the air. 
However, in the solar corona a much more suitable process is a source-region 
expansion which acts as the three-dimensional (3D) piston. If the expansion is 
impulsive enough it creates a large-amplitude perturbation, whose leading edge 
steepens due to non-linear effects, i.e. wave elements of higher amplitude move 
faster. Eventually, a discontinuity occurs in the wavefront profile, meaning that 
the shock is formed. Whereas the ID MHD piston problem (planar wave) can 
be solved analytically (Mann, 1995; Vrsnak and Lulic, 2000a), an analogous 2D 
or 3D problem can be treated analytically only by applying severe assumptions 
and approximations, and after all, a numerical evaluation is needed in any case 
(see, e.g., Zic et at, 2008; Afanasyev and Uralov, 2011). 

Thus, numerical MHD simulations are required to study a 2D and 3D piston 
mechanism for the magnctosonic-wave generation. Bearing in mind that the wave 
formation and evolution are strongly influenced, or even dominated, by physical 
properties of the environment and the characteristics of the driver itself, there 
are two alternatives in approaching this complex physical problem. One way is to 
set up the initial conditions as closely as possible to the real situation in which a 
particular wave has occurred, and to perform a full 3D simulation that provides a 
detailed quantitative analysis of the specific event. Such an approach, providing 
detailed insight into the physics behind a particular event, including coronal 
diagnostics, were performed by, e.g., Uchida, Altschuler, and Newkirk (1973), 
Wang (2000), Wu et al. (2001), Ofman and Thompson (2002), Ofman (2007), 
Cohen et al. (2009), Schmidt and Ofman (2010), Downs et al. (2011), Selwa, 
Pocdts, and DeVore (2012). Another way is to start from a somewhat simplified 
initial situation, which provides more extensive parametric studies and gives a 
more general view on the problem. In this type of simulation the CME is usually 
represented by an erupting 2D structure, anchored in the inert photosphere 
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(see, e.g., Chen et al., 2002, Chen, Fang, and Shibata, 2005, Pomocll, Vainio, 
and Kissmann, 2008, Wang, Shen, and Lin, 2009). 

In this article we consider some simple initial configurations/geometries to get 
an insight into the most basic characteristics of nonlinear processes governing the 
MHD wave formation and evolution in general. The idea is to isolate the basic 
processes that stand behind the wave formation in an idealized surrounding, 

1. e. to identify effects that are present regardless of the specific properties of 
the environment. In the follow-up article, more realistic configurations will be 
considered, including the chromosphere/corona density and Alfven speed profile, 
magnetic field line-tying, and the arcade expansion accompanied by an upward 
motion. These more advanced simulations will be compared with the results 
presented in this article, which will help us to distinguish the effects that are 
intrinsic to the MHD wave formation from those governed by the environment. 

Special attention is paid to the wavefront steepening, i.e., the shock formation 
process, in a planar and cylindrical geometry. For the simulations we employ the 
Versatile Advection Code (VAC: Toth, 1996; Goedbloed, Keppens, and Poedts, 
2003). This numerical code was developed at the Astronomical Institute at 
Utrecht, in collaboration with the FOM Institute for Plasma Physics, the Math- 
ematics Department at Utrecht, and the Centrum Wiskundc and Informatica 
(CWI) at Amsterdam. It is suitable for the analysis of a broad spectrum of 
astrophysical phenomena, including magnctohydrodynamic (MHD) shock waves. 

In Section 2 the model employed and the simulation procedure are briefly 
described. In Section 3 we present the results, first considering a planar geometry, 
so that the outcome can be compared with the analytical results, and then 
switching to a cylindrical geometry, which is more closely related to a coronal- 
arcade eruption or a coronal-loop expansion. In Section 4 we discuss the results 
and compare them with observations. 

2. The Model 

In the following, we consider perpendicular magnctosonic waves, where we focus 
on a planar and cylindrical geometry. This allows us to set the magnetic field in 
the z-direction, whereas the x and y magnetic-field components, as well as the 
z-component of the velocity, are always kept zero (B x = 0, B y = 0. v z = 0). 
Furthermore, all quantities arc invariant along the z-coordinate, i.e. we perform 
2.5D simulations, where the input and the basic output quantities are the density 
[p] the momentum [m x = pv x , m y = pv y ) and the magnetic field [B z ]. Note 
that although we perform 2.5D simulations, physically it is a one-dimensional 
problem. 

We use a two-dimensional [2D] numerical mesh containing 995 x 995 cells, 
supplemented by two ghost-cell layers at each boundary, which are used to 
regulate the boundary conditions (thus, the complete grid consists of 999 x 999 
cells). We apply continuous boundary conditions, meaning that gradients of all 
quantities are kept zero by copying the variable values from the edge of the 
mesh into the ghost cells. All quantities are normalized, so that distances are 
expressed in units of the numerical-box length [L = 1], velocities are normalized 
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to the Alfven speed [va] , and time is expressed in terms of the Alfven travel time 
over the numerical-box length \pA = L/va]- We apply the approximation j3 = 0, 
where /3 is the plasma-to-magnctic pressure ratio. The origin of the coordinate 
system is set at the numerical-box center. 

We will consider two basic initial configurations, resulting in a planar wave 
and a cylindrical one. In the planar option, all quantities are invariant in the 
y-dircction, i.e. quantities depend only on the x-coordinatc. In the cylindrical 
option, all quantities depend only on the r-coordinate, where r 2 = x 2 + y 2 . 

In all runs, we set up the simulation as an initial-value problem, starting 
from an unstable magnetic-field configuration of the source region. Thus, the 
source- region expansion is not fully under control, i.e. we do not prescribe the 
time-profile of the "driver" motion. The overall characteristics of the source- 
region expansion (the acceleration impulsiveness and the maximum speed) are 
regulated only indirectly, by increasing or decreasing the initial force imbalance. 
More precisely, we start from an initial configuration where the force balance is 
distressed by the excess magnetic pressure, i.e. the space-time evolution of the 
plasma flow is entirely determined by the initial spatial profile of B 2 . 

Specifically, we set a "parabolic" profile for the initial magnetic field within 
the source region: 



where Bo represents the magnetic field at x = and b(x) defines the field strength 
profile within the source region. We employ the form b = (Bq — B 2 )/xq, where 
xq is the initial source-region size and B e represents the external magnetic field 
strength outside the source region. In the cylindrical configuration we use the 
same function, only replacing x by r. The initial magnetic-field profile is drawn 
in Figure la by a red line. For the initial source-region size we take xq = 0.1; 
beyond x = xq we set B e = 1 and p e = 1. To make the source region more inert, 
and to better visualize the source region, we increased the density within the 
source region to p = 2 (see the red line in Figure lc). 

At the beginning, the plasma is at rest, v — (see the red line in Figure le). 
The considered profile of B z is characterized by a magnetic-pressure gradient 
[d(B 2 /2fio)/dx] which causes the initial outward acceleration. The acceleration 
increases linearly from at x = to a maximum value at the source-region 
boundary (hereinafter denoted also as "contact surface" , or "piston" ) . The mo- 
tion of the source-region boundary is tracked by identifying a surface within 
which the mass content equals the initial one. 

In the cylindrical geometry we also use another initial magnetic-field profile 
for the source region. It has the form: 



where ro represents the source-region size, and B e is the magnetic-field strength 
for r > ro . In this case the magnetic pressure gradient and the initial acceleration 
are zero at the source-region center and at the source surface, whereas the peak 




(1) 
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value is attained within the source-region body, at r/ro = 7r/4. The density is 
again set to p = 2 within r < ro , and p = 1 for r > r$ . 

Although we do not intend to reproduce directly any specific coronal struc- 
ture, note that the hereafter presented evolution of the two cylindrical configura- 
tions employed can depict, to a certain degree, the coronal- wave formation caused 
by a lateral expansion of the source region placed in a vertical magnetic field of 
a quiet corona. Such an expansion can occur, e.g. in the impulsive-acceleration 
stage of CMEs (the so-called "lateral overexpansion" ; see, e.g., Kozarcv et al., 
2011; Patsourakos, Vourlidas, and Stenborg, 2010), or presumably, in legs of 
impulsively heated flaring loops. 

3. Results 

3.1. Formation and Propagation of a Planar Shock 

First we analyze the formation and propagation of a perpendicular shock in the 
planar geometry. The aim is to compare the numerical results with the analytic 
theory for planar MHD waves (Vrsnak and Lulic, 2000a; 2000b) and to have 
the reference-results when the influence of the geometry on the results will be 
considered. 

The formation and propagation of the wave is presented in Figure 1 . In Fig- 
ures la and b we show the magnetic-field profiles [B z {x)\ in Figures lc and d the 
density profiles [p(x)] and in Figures le and f the flow-speed profiles [v x (x)}. The 
graphs in the left column show the formation phase of the wave, whereas those 
on the right side represent the propagation phase. The kinematics of various 
features recognized in the density profiles in Figures lc and d (the wave front 
leading edge, the wave peak, a dip between the wave and the piston, and the 
source- region boundary) are displayed in Figure 2. 

Due to the magnetic pressure gradient of the unstable initial configuration, 
the source-region expansion starts immediately at t = 0. The acceleration is 
strongest at the source-region surface, whereas the source-region center (x = 0) 
stays at rest. Over most of the source-region volume, the density decreases due to 
the expansion, whereas close to the contact surface it increases due to the velocity 
gradient. The kinematics of this density peak closely follows the kinematics of 
the contact surface, just slightly lagging behind it. A peak density, p = 2.34, is 
attained around t = 0.06. The source-region boundary accelerates until t s=s 0.1, 
attaining a speed of v = 0.4. After that, it gradually slows down, and stops 
around t w 0.35 (see kinematics shown in Figure 2a). During the accelerated- 
expansion phase, the flow speed increases, attaining a value of v ~ 0.4 around 
i w 0.1 (Figure le), i.e. the fastest flow elements are adjusted to the piston 
motion. 

Ahead of the contact surface, a wave front forms as a result of the source-region 
expansion. It can be easily recognized in the magnetic-field and density profiles 
shown in Figures la and c. The wave detaches from the source region after t ps 0.1 
(i.e. after the piston acceleration-phase ends), and continues to evolve as a freely 
propagating simple wave (for a hydrodynamic analog see Sections 101 and 102 
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Figure 1. Formation and propagation of a perpendicular shock in the planar geometry: spatial 
profiles of the magnetic field (a, b), density (c, d), and flow speed (e, f). Left panels (a, c, e): 
the beginning of the wave formation (increasing wave amplitude); right panels (b, d, f): the 
shock formation phase (steepening of the wavefront profile). The initial magnetic field in the 
source-region center is Bo = 2. Times are displayed in the inset. All quantities presented are 
normalized: distance x is expressed in units of the numerical-box length [L = 1], velocity v x 
is normalized to the Alfven speed [va], and time t is expressed in terms of the Alfven travel 
time over the numerical-box length [t^ = L/va]. 
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Figure 2. Kinematics of various wave features and the source-region boundary (thin-solid 
line — the wavefront leading edge; thick-solid line - the wave crest; dot-dashed line - a trailing 
density dip; dotted line - the source-region boundary): a) distance versus time; b) velocity 
versus time. 



in Landau and Lifshitz, 1987). Note that a dip in the density profile, formed 
between the wave peak and the contact-surface peak, never gets values p < 1. 
On the other hand, the density in central parts of the source region becomes 
strongly depleted. 

The wavefront steepens in time, whereas its amplitude remains constant, 
staying at values of p = 1.44, B = 1.44, and v = 0.42, respectively. The shock 
formation is completed at t re 0.26. 

The kinematics of the wave leading edge, the wave peak, a rarefaction dip, and 
the piston, measured from the density profiles shown in Figure 1, are displayed 
in Figure 2, revealing that the piston accelerates until t re 0.08. Thereafter, it 
continues to move at an approximately constant speed of v re 0.4 until t re 0.13. 
During this period the wave amplitude increases (see Figure 1) and the wave- 
crest phase speed increases from w re 1 to w — 1.76. At the same time, the 
wavefront leading edge moves at w re 1. The wave crest reaches the leading edge, 
i.e. the shock formation is completed, around t re 0.25. After that, the shock 
front moves at a speed of w = 1.35, consistent with the Rankine-Hugoniot jump 
relations. The described evolution of the source/wave system and its kinematics 
is fully consistent with the analytical model presented by Vrsnak and Lulic 
(2000a). 

After t re 0.13, the source-region expansion gradually decelerates, and prac- 
tically stops at t re 0.35. A density dip between the wave peak and the piston, 
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Figure 3. Relationship between the wave speed [w] and the flow speed [v] for the planar wave. 
Numerical results (red squares) are compared with the analytical relationship w = 1 + 3v/2 
(blue line) derived by Vrsnak and Lulic (2000). 

which forms around t ~ 0.11, closely follows the kinematics of the source-region 
boundary, being only slightly faster then the piston. 

We repeated the simulations using various values of Bq, to analyze how 
the evolution of the piston/wave system depends on the impulsiveness of the 
source-region expansion. A stronger Bq results in a more impulsive source- 
region acceleration, which leads to a higher shock amplitude and Mach number. 
Furthermore, the shock is created earlier and closer to the piston, so in the case 
of extremely impulsive expansions, the shock-sheath region and the source region 
cannot be clearly resolved. 

In Figure 3 we show the dependence of the phase speed [w] of the perturbation- 
segment at the wave crest (before being shocked) as a function of the correspond- 
ing flow speed v. In the graph we display the results for Bq = 1.5, 2.0, 3.0, and 
5.0. For Bq = 2 we also measured w and v at several suitable wavefront-segments 
ahead of the wave crest (the lowest v- values in Figure 3). The results are fully 
consistent with the outcome of the analytical theory for (3 — presented by 
Vrsnak and Lulic (2000a), where the relationship w = 1 + 3v/2 was established. 

Whereas in the non-shocked phase the wave behavior is consistent with the 
analytical theory, in the shocked phase for Bq > 3, corresponding to w > 1.8, 
the numerical results start to deviate from the analytical Rankine-Hugoniot re- 
lations, the disagreement increasing with the increasing amplitude. The equation 
of continuity and the relationship between Mach number and the downstream 
flow speed behave as expected, but the relationship between Mach number and 
the downstream/upstream density jump deviates from the analytical results. 
This is probably due to numerical effects and the fact that at very high values of 
Bq it becomes impossible to clearly resolve the compression at the source-region 
surface and the shock itself. 

3.2. Cylindrical Geometry 

3.2.1. Wave Formation 

In Figure 4 the formation and propagation of the wave in the cylindrical geometry 
is presented. Spatial profiles of the magnetic field [B z (r)\ are shown in Figures 4a 
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Figure 4. Formation and propagation of the perpendicular shock in the cylindrical geometry 
for the initial magnetic field profile given by Equation (1) with Bo = 2: spatial profiles of the 
magnetic field (a, b), the density (c, d), and the flow speed (e, f). Left panels (a, c, and e) show 
the beginning of the wave formation (increasing wave amplitude); right panels (b, d, and f) 
present the shock formation phase (steepening of the wavefront profile). Times are displayed 
in the inset. All quantities presented are normalized: distance r is expressed in units of the 
numerical-box length [L = 1], velocity v r is normalized to the Alfven speed [va], and time t is 
expressed in terms of the Alfven travel time over the numerical-box length [fyt = L/va]- 
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and b, the density profiles [p(r)] are presented in Figures 4c and d, whereas in 
Figures 4e and f the flow speed [iv(r)] is displayed. The initial magnetic- field and 
density profiles are defined in the same way as in the planar case (Equation (1)), 
only replacing x — > r. In the left column the wave formation phase is shown, 
whereas the right column represents the propagation phase. In Figure 5, the 
kinematics of various features recognized in Figure 4 are shown. 

As in the planar case, the source- region expansion starts immediately at t = 0, 
the acceleration being strongest at the source-region surface. The source-region 
center (r = 0) remains at rest at all times (Figures 4e and f). The density 
within the source region starts to decrease due to the expansion, whereas at the 
contact surface it increases due to the flow-speed gradient. The source-region 
expansion initially accelerates, attaining v = 0.28 at t « 0.07 (see Figure 5). 
Note that the acceleration phase is shorter than in the planar case, the peak 
speed is considerably lower, and the w« const, phase is absent. After attaining 
the maximum speed, the piston gradually slows down, stops around t ps 0.2, 
and then retreats slowly towards the initial position (see the kinematics shown 
in Figure 5). During the accelerated-expansion phase, the flow speed increases, 
reaching v — 0.28 around t ~ 0.04 (Figure 4e). Note that, unlike in the planar 
geometry, the plasma flow is not fully synchronized with the piston motion. 

Ahead of the contact surface, the wavefront forms as a result of the source- 
region expansion. It can be readily recognized in the magnetic field and density 
profiles shown in Figures 4a and c. The wave detaches from the source region 
around t w 0.08, having an amplitude in p and B of around 1.22. After that, 
the perturbation continues to propagate as a freely propagating simple wave, but 
unlike in the planar case, the amplitude of the wave decreases with distance (Fig- 
ure 4). The wavefront steepens with time, whereas the peak flow-speed decreases. 
A discontinuity in the leading-edge profile occurs, i.e. the shock formation begins 
at t ps 0.15. The shock is fully completed at t pa 0.28, when it has an amplitude 
of p = 1.16 and v = 0.16. 

Note that flows within the source region are more complex than in the planar 
configuration. We also stress that a dip in the density profile, which forms 
between the wave peak and the contact-surface, now deepens to a value of 
p = 0.88, i.e. the rarefaction region forms (p < 1), as in the case of cylindrical 
hydrodynamic waves (see Section 102 in Landau and Lifshitz, 1987 and references 
therein). 

3.2.2. Wave Kinematics 

The kinematics of the piston and the wave, estimated from the density profiles 
displayed in Figure 4, are shown in Figure 5. Comparing Figures 4 and 5 one 
finds that during the piston acceleration the wave amplitude first increases, but 
then starts to decrease, even before the piston reaches its maximum velocity. 
The phase speed of the wave crest increases from w w 1 to w ~ 1.3 at f > 0.15, 
thereafter it gradually decreases. 

A dip between the wave peak and the piston, which forms around t ps 0.08, 
first closely follows the piston kinematics, but then, after the dip becomes char- 
acterized by p < 1 at t pa 0.2, it "detaches" from the piston and attains a speed 
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Figure 5. Kinematics of the wave and the source-region boundary derived from the density 
profiles shown in Figure 4: a) distance versus time; b) velocity versus time. Thin-solid line — the 
wavefront leading edge; thick-solid line — the wave peak; dashed line — a deep front measured 
at p = 1; dot-dashed line - the dip minimum; dotted line — the source-region boundary. All 
quantities presented are normalized: radial distance r is expressed in units of the numerical-box 
length [L = 1], velocities v and w are normalized to the Alfven speed [va], and time t is 
expressed in terms of the Alfven travel time over the numerical-box length [tj\ = L/va\- 



of w « 1 around t « 0.35. Note that a segment of the dip characterized by p = 1 
moves at a speed of w w 1 all the time. 

The relationship between the shock speed and the downstream flow speed 
in the cylindrical geometry is more complex than in the planar case. This is 
illustrated in Figure 6a, where the downstream flow-speed [v] is shown versus 
the phase-speed of the wave crest [w] . Analogously, in Figure 6b we present the 
dependence of the downstream peak density [p] on the phase speed [w]. Note 
that the displayed values are based on smoothed curves w(t), v(t), and p(t). 
The presented graphs show that initially, during the wave formation phase, the 
wave phase-speed increases, while the amplitudes of v and p are almost constant, 
showing only a slight increase. Then, the wave speed [w] remains almost constant, 
whereas the wave amplitude decreases. Eventually, in the third step, both the 
wave propagation speed and its amplitude decrease. 

The highest values of v and p are attained roughly at a time when the 
maximum speed of the piston is reached. The shock formation starts, i.e. a 
discontinuity occurs at the leading edge of the wave profile, around the "nose" 
of the curves presented in Figure 6, which also approximately coincides with the 
end of the piston expansion. Thus, roughly speaking, the upper branch of the 
v(w) and p(w) curves corresponds to the "driven phase" of the wave, whereas 
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Figure 6. The evolution of: a) downstream flow-speed amplitude versus phase-speed of 
the wave crest; b) density amplitude versus phase-speed of the wave crest. Solid-blue and 
dashed-red lines show results for the initial configuration defined by Equations (1) and (2), 
respectively, with Bo = 2. Arrows indicate the course of the temporal sequence. Velocities v 
and w are normalized to the Alfven speed [v^]. 

the lower branch represents the "decay" of a freely-propagating simple wave. 
The shock is completed when the wave speed attains a value ofio« 1.11 at the 
lower branch of the curve. 

3.2.3. Piston Impulsiveness 

We repeated the procedure presented in Sections 3.2.1 and 3.2.2 for several 
values of the maximum magnetic-field strength in the source-region center [Bq] 
to inspect the role of the impulsiveness of the piston expansion. In particu- 
lar, we applied Bq = 1.1, 1.5, 2.0, 3.0, and 5.0. A stronger field Bo causes a 
more impulsive acceleration of the piston, which results also in a higher wave 
amplitude and wave-crest speed, and consequently, an earlier formation of the 
shock. On the other hand, the evolution of the system, as well as the relationship 
between different parameters, docs not depend qualitatively on the impulsiveness 
of the piston acceleration. Morphologically, the main difference between very 
impulsive piston accelerations and more gradual ones is that in the former case 
the shock forms very close to the contact surface. Because of this, in the case of 
a very impulsive source-region expansion, it is not possible to follow the shock- 
formation phase, since the wavefront and the piston cannot be resolved. On the 
other hand, we note that for Bq = 1.1 and 1.5 the shock did not form within 
the numerical box, which implies that in reality, particularly bearing in mind 
dissipativc effects, the coronal shock would not be formed if the source-region 
acceleration is not impulsive enough. 
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Figure 7. a) Peak phase-speed [w] of the wavefront (solid line and diamond symbols) and 
the speed [v] of the piston (dotted line and crosses) for five values of Bo- b) The phase-speed of 
the wavefront around the peak amplitude (squares, triangles, crosses, asterisks, and diamonds; 
the values of Bo are given in the inset) presented as a function of corresponding flow speed. 
The solid line represents the peak velocity of the wave as a function of the maximum speed 
of the piston for the same five values of Bo . Velocities v and w are normalized to the Alfvcn 
speed [va]. 

In Figure 7 a the peak velocity of the wavefront is compared with the peak 
velocity of the piston for different values of Bq. In the considered range, the wave 
speed is much larger than the piston speed, so the distance from the wavefront 
and the piston rapidly increases. However, the graph shows that beyond Bq « 1.5 
maximal piston velocities are proportional to Bq, whereas the wave speed shows 
a nonlinear trend, i.e. the slope of the w(Bq) curve gradually decreases. This 
implies that for a very impulsive expansion of the source region one can expect 
that the velocities of the shock and the piston become comparable, and that the 
separation is small. 

In Figure 7b we present the dependence of the maximum speed of the wave 
crest on the maximum speed of the corresponding downstream flow speed for 
all considered values of Bq. The displayed data points arc numerical values from 
the "nose" of the w(v) curves (non-smoothed), analogous to the one shown in 
Figure 6a for Bq = 2. Peak values of w, estimated from smoothed w(v) curves, 
are presented as a function of the maximum speed of the piston by a solid line. 
Note that this "piston-curve" is shifted to the right with respect to the presented 
data points, implying that the piston speed is somewhat higher than the flow 
speed. The difference increases with the increasing piston speed, i.e. with the 
impulsiveness of the expansion. 
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The main feature of Figure 7b is that the relationship between the wave speed 
and the downstream flow speed [w = l+3v/2] is not valid in cylindrical geometry. 
The w(v) dependence is not linear, but is closer to a power-law. A least-square 
fit of the form w — 1 = av b gives w = 1 + 0.9 v 0A5 , with a correlation coefficient 
of R = 0.91. On the other hand, the relationship between the maximum wave 
speed and the maximum piston speed is well described (R = 0.99) by the function 
w = 1 + 1.26 v 1 ^. 

3.2.4-. Initial Configuration 

To check how much the initial magnetic-field structure in the source region 
affects the process of the wave formation and evolution, we also applied the 
magnetic-field configuration defined by Equation (2). In this case, the steepest 
magnetic-pressure gradient is located within the source region, i.e. not at its 
edge as was the case described by Equation (1). The outcome of the simulation 
for Bo = 2 is presented in Figure 8. 

The analysis of the data in Figure 8 shows that there are no significant 
differences in the overall wave kinematics for the two configurations considered 
(thus, the corresponding graphs are not shown). To illustrate the similarity of 
the two kinematics, we also included in Figure 6 the results concerning the wave 
formation/propagation resulting from the configuration defined by Equation (2). 

However, Figure 8 reveals a considerable difference in the morphology of 
the evolving piston/wave system. The strongest magnetic pressure gradient is 
initially located within the source region, whereas it is zero at its edge. Thus, the 
highest acceleration occurs at r < ro, and consequently, the initial compression 
forms within the source region. This causes a more complex flow pattern within 
the source and obscures the initiation of the wave, i. e. the source- region bound- 
ary and the wave cannot be clearly distinguished, and the term piston becomes 
meaningless. 

The wave leading edge leaves the source region at t — 0.015, whereas the 
wave crest detaches from the source-region boundary at t = 0.08, when it has 
the highest amplitude (p = 1.22). The shock formation starts at t « 0.19 and 
is completed at t « 0.32, when it has an amplitude of p = 1.13 and v = 0.13, 
i.e. the formation is delayed, and the amplitude is lower compared to the case 
based on Equation (1). Similarly, the dip between the wavefront and the source 
region is somewhat shallower (p = 0.92). 

4. Discussion and Conclusion 

In this article we presented numerical simulations of the formation and evolution 
of large-amplitude MHD simple waves. We considered the very basic initial 
configurations to educe general characteristics of the MHD shock formation in 
an idealized homogeneous environment. The main purpose is to have reference 
results that can be compared to the results of more sophisticated simulations 
that consider more realistic characteristics of the environment. Now that we 
have at our disposal the results of the simulations presented in this article, such 
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Figure 8. Formation and propagation of the perpendicular shock in the cylindrical geometry 
for the initial magnetic-field profile given by Equation (2), with Bo = 2: spatial profiles of 
the magnetic field (a. b). density (c, d), and flow speed (e, f). Left panels (a, c, and e) show 
the beginning of the wave formation (increasing wave amplitude); right panels (b, d, and f) 
present the shock formation phase (steepening of the wavefront profile). Normalized times are 
displayed in the inset. 
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a comparison will help us to distinguish which characteristics are a consequence 
of basic processes and which are caused by the details of a given environment. 
Furthermore, it should be noted that in spite of the high level of idealization, 
the cylindrical configurations employed can represent, to a certain degree, the 
initial stage of the coronal-wave formation by the lateral expansion of a CME 
during its impulsive-acceleration stage, or presumably, by the expanding leg of 
the flaring loop. 

The most general outcome, common to all situations analyzed, is that a more 
impulsive source- region expansion results in a shorter time/distance needed for 
the shock formation, consistent with analytical considerations (e.g. Vrsnak and 
Lulic, 2000a; Vrsnak and Lulic, 2000b; Vrsnak, 2001; Zic et al, 2008) and obser- 
vations (Vrsnak, 2001). The simulations show that in the most impulsive events 
a shock forms very close to the source-region boundary and it is initially difficult 
to resolve the two entities. This explains why in some studies the coronal EUV 
waves are (erroneously) identified as CME flanks (for a discussion see Cheng 
et al, 2012). On the other hand, when the piston acceleration is low, the wave 
amplitude remains small and the wavefront steepening is slow. Thus, weakly 
accelerated eruptions are not likely to result in an observable coronal wave. 

For the case of a planar magnetosonic wave we have confirmed the relationship 
between the wave speed and the flow speed [w — 1 + 3v /2] that was analytically 
derived by Vrsnak and Lulic (2000a). At small amplitudes, the numerical simula- 
tions reproduce the Rankine-Hugoniot jump relations after the shock formation 
is completed. However, at large amplitudes the numerical results deviate from 
the analytical theory, most likely due to the numerical resolution. 

From the observational point of view, the cylindrical geometry is far more 
interesting, since it can give us insight into the process of shock formation 
caused by a magnetic-arcade expansion, including the amplitude fall-off due to 
energy conservation (Zic et al., 2008). In this article we have analyzed only the 
most general characteristics of the perpendicular-shock formation, generated by 
a flux-tube expansion in an idealized homogeneous environment. Such a process 
represents a two-dimensional piston mechanism of the shock-wave formation. 
The basic difference from the planar case (one-dimensional piston) lies in the fact 
that in the cylindrical geometry there arc two competing effects involved in the 
shock-formation process. One is the nonlinear steepening of the wavefront profile 
(as in the planar geometry) , and the other is a decrease of the wave amplitude 
with distance due to energy conservation (absent in planar geometry). 

Of course, spherical geometry would be more relevant than the cylindrical 
option, since real EUV waves expand spherically. This would certainly modify 
the results, since the decrease rate of the wave amplitude would be governed by 
the oc r~ 2 effect rather than oc r _1 . This aspect was treated in detail by Zic et al. 
(2008) in a semi-analytical study which showed that the shock formation time 
and distance depend much more on the characteristics of the piston acceleration 
than on the choice of the geometry. The geometrical effect (i.e. the oc r~ 2 aspect 
of the energy conservation) becomes dominant only after the acceleration phase. 
Note that the same conclusion can be drawn from the results presented in this 
article by comparing the outcome for the planar and cylindrical case. In this 
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respect, let us note that Zic et al. (2008) have also shown that the value of the 
plasma-to-magnetic pressure ratio [/3] does not-play a significant role as well. 

In our study we considered two different types of the initial configuration: one 
resulting in the highest initial acceleration at the source-region boundary, and 
another, causing the strongest acceleration within the source-region body. The 
performed simulations show that, although there are differences in the evolution 
of the source region, the process of the shock-wave formation docs not differ 
significantly, i. e. the evolution of the perturbation and the wave kinematics are 
similar in both cases. 

The most important outcome of the analysis is that the formation of a perpen- 
dicular MHD shock is expected already for relatively low expansion velocities, 
as low as 10-20% of the Alfven speed. This implies that a lateral expansion of 
the eruption in the early phase of CMEs is a viable mechanism of the coronal 
wave formation (Kienreich, Temmcr, and Vcronig, 2009; Patsourakos, Vourlidas, 
and Stcnborg, 2010; Muhr et al, 2010; Vcronig et al, 2010; Grechncv et al., 
2011; Kozarcv et al, 2011; Liu et al, 2012; Temmcr, Vrsnak, and Vcronig, 2012). 
Furthermore, our simulations show that at the beginning of the wave formation 
it is difficult to distinguish the wave and the source-region expansion, especially 
for the case defined by Equation (2), where the strongest acceleration occurs 
within the source region. 

The presented analysis shows that the wave initially accelerates from w > va 
to a maximum phase speed, which depends on the impulsiveness of the source- 
region expansion. In the decay phase, the wave-crest velocity decreases, w — > va- 
Thus, the initial and the late phase of the coronal wave could be used for the 
coronal diagnostics, since measurements of the wave kinematics in the acceler- 
ation and deceleration phase should reflect the coronal Alfven speed. Similarly, 
the traveling density depletion that forms in the wake of the wave travels at 
w ~ va- Such depletions are sometimes observed in the base-difference or base- 
ratio EUV images, appearing as traveling coronal dimming behind the wavefront 
(Thompson et al, 2000; Chen et al, 2002; Zhukov and Auchcre, 2004; Muhr 
et al, 2011). Thus, such features can be also used to estimate va in the quiet 
corona. Finally, we note that after the acceleration stage, the compression region 
associated with the source-region boundary becomes a stationary feature. This 
might be related to stationary brightenings that are sometimes observed behind 
the outgoing wave (Muhr et al, 2011). 

In cylindrical geometry, the wave amplitude and the wave phase-speed are 
related in a relatively complex manner. In the "driven phase" , the wave ampli- 
tude at a given wave speed is higher than in the "decay phase" . In the transition 
between these two phases, the phase speed is almost constant for a certain period 
of time, while the wave amplitude decreases. This results in a loop form of 
the p(w) and v(w) evolutionary curves, which is consistent with observations 
presented by Muhr et al (2013), where the dependence of the wave amplitude 
on the wave speed forms a closed, hysteresis-like, curve. 

In Section 3 we have presented results in a normalized form, where velocities 
are expressed in units of the Alfven speed and the time is expressed in units 
of the Alfven travel time across the numerical box. As an illustration, let us 
assume that the diameter of the coronal source region is 2ro = 100 Mm, which 
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implies that the numerical box corresponds to L = 500 Mm, since we used 
ro = 0.1. If we assume that the Alfven speed in the quiet corona is in the range 
of da = 250 kms -1 (Warmuth and Mann, 2005), we obtain for the Alfven travel 
time tA = L/va = 2000 seconds. The same would be obtained for, e.g. L = 1000 
Mm and i>a = 500 kms" 1 . 

Applying these values, one finds that the wave typically forms and steepens 
into a shock a few minutes after the onset of the source-region expansion. The 
delay is shorter for a higher source-expansion velocity, i.e. for a higher wave 
velocity. The distance at which the wave crest forms and detaches from the 
source- region boundary, i. e. the distance at which the wave should become ob- 
servable, is in the range rj 100-200 Mm. Such time delays and starting distances 
are fully consistent with observations of Moreton waves (e.g. Warmuth et al., 
2004b), type II bursts (e.g. Vrsnak et al, 1995), and EUV waves (e.g. Kienreich 
et al, 2011; Liu et at, 2012). 

To conclude, the simulations presented show that already the simplest initial 
configurations are able to reproduce most of the features observed in typical 
large-amplitude, large-scale coronal waves, including the morphology, kinemat- 
ics, and scalings. Our next step will be to perform similar simulations, but 
employing more realistic initial configurations that depict a magnetic arcade 
anchored in the photosphere, and include a realistic density profile of the chro- 
mosphere, transition-region, and corona. This will enable us to distinguish the 
effects that are intrinsic to the MHD wave formation from those governed by 
the environment. 
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